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Abstract. We prove that planar homeomorphisms can be approxi- 
mated by diffeomorphisms in the Sobolev space W^'"^ and in the Royden 
algebra. As an application, we show that every discrete and open planar 
mapping with a holomorphic Hopf differential is harmonic. 



1. Introduction 

It is a fundamental property of Sobolev spaces y^^''*', l^p<oo, that any 
element can be approximated strongly (i.e., in the norm) by '^'^ smooth 
functions, or by piecewise affine ones. In the context of vector- valued 
Sobolev functions, that is, mappings in W^^'^i^^ M"), invertibility comes into 
play. Indeed, the studies of invertible Sobolev mappings are of great impor- 
tance in nonlinear elasticity [21 [T^ [26l [35] . The following natural question 
was put forward by John M. Ball. 

Question 1.1. [3] If u G #^^''^($1, M") is invertible, can u be approximated 
in ^^'P by piecewise affine invertible mappings? 

J. Ball attributes this question to L.C. Evans, who was led to it through 
his investigation of the partial regularity of minimizers yL3j of neohookean 
energy functionals O O [TJ [33] . We provide an affirmative solution of the 
Ball- Evans problem in the case p = n = 2. The most general formulation 
of our result administers Royden algebras and £/o{^), see Section [2j 

We write 

£[h] = £n[h] := \\Dh\\%,^^^ = JjDh{z)\' dz 
where \Dh\ is the Hilbert- Schmidt norm of the differential. 

Theorem 1.2 (Approximation by diffeomorphisms). Let h: Q be a 

homeomorphism of Sobolev class /^^^ (0,0*). Then for every e > there 
exist a diffeomorphism H : Q, Q,* such that 
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(i) H -he £/o{^) 
(ii) \\H - < e 

(Hi) £[H] i^£[h]. 

Part (iii) is nontirivial only in the finite energy case, <Sq[/i] < oo. Let us 
note that the existence of smooth approximation implies the existence of 
piecewise-affine approximation, since a diffeomorphism can be triangulated. 
(In the converse direction, a piecewise-affine mapping can be smoothed in 
dimensions less than four [25 1, but not in general.) Partial results toward the 
Ball-Evans problem were obtained in [23] (for planar bi-Sobolev mappings 
that are smooth outside of a finite set) and in [6] (for planar bi- Holder 
mappings, with approximation in the Holder norm). The articles [H [32] 
illustrate the difficulty of preserving invertibility in the process of smoothing 
a Sobolev homeomorphism. 

We also give an application of Theorem 11.21 to a problem that originated 
in a series of papers by Eells, Lemaire and Sealey [11[ [T2] 131] . It concerns 
the nonlinear differential equation 

(1.1) ^.{h-J^)=Q 

for mappings defined in a domain in the complex plane C. Naturally, the 
Sobolev space W^f{Vl^C) should be considered as the domain of definition 
of equation (II. ip . This places hzhz in Jfj^^(r2), so the complex Cauchy- 
Riemann derivative ^ applies in the sense of distribution. By Weyl's lemma 
hzhz is a holomorphic function. 

The expression Qh '■= hzhzdz® dz is known as the Hopf differential of h 
(named after H. Hopf, who employed a similar device, e.g., in [T9l Chapter 
VI] ) . It is clear that Qh is a holomorphic quadratic differential whenever h is 
harmonic, which is a general fact about energy-stationary mappings between 
Riemannian manifolds [10, (10.5)], [2T] and |34j. Eells and Lemaire inquired 
about the possibility of a converse result, e.g., for mappings with finite 
energy and almost-everywhere positive Jacobian |11[ (2.6)]. In this setting 
a counterexample was provided by Jost [20], who also proved the existence 
of ^'I'^-solutions of ([HJ in every homotopy class of mappings between 
compact Riemann surfaces. A more restricted form of the Eells-Lemaire 
problem, [12l (5.11)] and [3l], imposed the additional assumption that h 
is a quasiconformal homeomorphism, and was settled by Helein |T7] in the 
affirmative. Here we dispose with the quasiconformality condition and treat 
general planar homeomorphisms of finite energy. Since the inverse of such a 
homeomorphism need not be in any Sobolev class [18], some difficulties are 
to be expected. They shall be overcome with the aid of our approximation 
theorem 11.21 

Theorem 1.3. Every continuous, discrete and open mapping h of Sobolev 
class yf^^ {^^C) that satisfies equation p.ip is harmonic. 
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The failure of Theorem 11.31 for uniform hmits of homeomorphisms should 
be mentioned. This is illustrated by Example 14.11 

2. Background 

Let be a bounded domain in M? ~ C, nonempty open connected set. We 
consider a class (0,) of uniformly continuous functions /i : i7 — )• C having 
finite Dirichlet energy, and furnish it with the norm 

l|/i|U(Q) = \\hy{n) + \\Dh\\^2^Q-j < oo 

is a commutative Banach algebra with the usual multiplication of 
functions in which ||/ii/i2||£/(r2) ^ ll^ilLc/(r2)ll^2||js'(c)- The closure of "^0°° ($7) 
in will be denoted by =2/0 (^^)- Suppose, to look at more specific sit- 

uation, that = U is a Jordan domain; that is, a simply connected open 
set whose boundary F = dV is a closed Jordan curve. By goodness of the 
Caratheodory extension theorem p7l p. 18], there is a homeomorphism 
if: B U of the closed unit disk D = {,^: |^| ^ 1} that is conformal in D. 
After the change of variable, z = ^p{^), we obtain a function H{^) = 
in ^(B). The operation 

so defined is an isometry; ||T^/i||^(p) = ||/i||^(u)- Furthermore, 

Proposition 2.1 (A generalization of Poisson's formula). Let V be a Jordan 
domain. There is (unique) bounded linear operator 

Pu: ^(U) ^ 

such that 

Pu-Id: =e/(U) ^^o(U) 
A o Pu = 

We name Pu the Poisson operator. The energy of Puj/i does not exceed 
that of h. This fact is known as Dirichlet' s principle 

The proof of this proposition reduces to the case when U = B, by conformal 
change of variables. A routine verification of this case is left to the reader. 
We only indicate that the less familiar property Pnj/i — h £ J2^(D), for 
h e (B), needs to be justified. 

Corollary 2.2 (Harmonic replacement). Let Vt be a domain in C and U C 
U C a Jordan domain. There exists (unique) bounded linear operator 



Ru: s^{^) 
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such that, for every h £ .s/ (Q) 

{Ru/i = h 
ARu/i = 



on 
in U 



The Laplace equation yields £n\J^vh] ^ <frj[/i]- Equality occurs if and only if 
h is harmonic in U. 

A short proof of this corollary runs somewhat as follows. The unique 
harmonic extension of h: dV — t- C inside U given by Pu^ has the property 
that PiLj/i — h € J2^^(U). Therefore, the zero extension of Pu^ — h outside U, 
denoted by [Pu^ — h]o, belongs to £^(17). We define 



The desired properties of the operator Ru so defined are automatically ful- 
filled in view of Proposition 12.11 

Proposition 2.3. Let i7 6e a domain in C and U C U C 1^ a Jordan 
domain. Suppose that h £ £/ (fi) is a homeomorphism of onto h{Q) and 
h{V) is convex. Then Ru/i is homeomorphism in Q and is a harmonic 
diffeomorphism in U. 

The injectivity of Ru^ is the content of the Rado-Kneser-Choquet The- 
orem [21 p. 29]. Furthermore, planar harmonic homeomorphisms are 
smooth diffeomorphisms according to Lewy's theorem [Oj p. 20]. 

3. Smoothing Sobolev homeomorphisms. Theorem 11.21 

Proof of Theorem I i.^L We may and do assume that h is not harmonic, since 
otherwise H = h satisfies the desired properties by Lewy's theorem (men- 
tioned above). Let Zq £ be a point such that h fails to be harmonic in 
any neighborhood of Zo- By choosing the origin of the coordinate system 
we ensure that h[zo) does not lie on the boundary of any dyadic squares 
associated with the coordinate system. 

Let us choose and fix any e > 0. The construction of H proceeds in 5 
steps. We construct homeomorphisms h^ : ^ r^*, /c = 0, . . . , 5 such that 
Hq = h, hk € hk-i + /15 is a diffeomorphism, and \\hk — is 

bounded by a multiple of e for each k. In each step we modify the previous 
construction to gain better regularity. In steps 1, 2 and 4 we use harmonic 
replacement according to Proposition 12.31 In steps 3 and 5 we smoothen the 
mapping near the boundaries of the domains in which harmonic replacement 
was performed. The result of each step is denoted by /ii, . . . , /15. The finite 
energy case h £ £/ (0) requires a few additional details, which are provided 
at the end of each step. 

We begin with a decomposition of the target domain 



Ru/i := [Pu^ - /i]o + /i e 



00 



(3.1) 




i/=i 
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into closed nonoverlapping dyadic squares Qu C ft*. This decomposition 
is made by selecting the maximal dyadic squares that lie in il*. Thus the 
cover of by such squares is locally finite. The preimage of Q;^ under h, 
denoted by Uj,, is a Jordan domain in 0,. Hereafter Vi, will be referred to as 
the curved-square. In fact to every partion of Q* into closed squares there 
will correspond a partition of Q into closed curved-squares via the mapping 
h: Q Jl*, for example: 



Step 1. For each U^, we replace h: U^, Qj, with a piecewise harmonic 
homeomorphism hi : Q,^ that coincides with h on d\J,y. To this effect 

we partition the square Qjy, 

(3.2) C = QiuQ2u---UQ^, (n = n^ = 4'='0 

into congruent dyadic squares Qj,, i = 1, . . . ,n. The number n, depending 
on v, will be determined later. For the moment fix and look at the 
homeomorphisms 

h:Vl = h-\Wu)°-^Wu 
These mappings belong to the Royden algebra £/(V^^). With the aid of 
Propositions 12 . 1 1 and 12 . 31 we replace each h: Ul ^ Ql with a harmonic home- 
omorphism hi: UJ, that coincides with h on dUl, i = 1,2, ... ,n. 
Such mappings are ^°°-smooth diffeomorphisms : QJ,. Moreover, 
hi-hG s^o{Ui) and 



(3.3) 



^ViK] ^ ^Viih] forl,2,...,n 

^au' [^v] = ^9U' [h] 5 because hi = h on dUI 



We obtain a piecewise harmonic homeomorphism by gluing hi together along 
the common boundaries of VI . Denote it by 



hi £h + £^oiVu) 



Precisely we define 



h1 = h + J2[K-h]o 

i=l 

Here and in the sequel the notation [ip]o for ip G =e/o(U) stands for zero 
extension of (p to the entire domain 17. Obviously [ip]o G £/oi^)- The above 
construction depends on the number n. For v fixed we actually have a 
sequence {/i"}ri=i,2,... that is bounded in J2/(U^). However, we have uniform 
bounds independent of n, 

WKWviU;;) ^ diamQ^ 
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and 



The key observation is that 



(3.4) 

Kn—^oo 

Indeed, for z G U* we have 



\K{z) - h{z)\ ^ diamQt = ^diamQ,, 

v ^ 

Thus h" ^ h uniformly on as n ^ oc. On the other hand the differential 
matrices Dh^ are bounded in ^^(Ujy, M'^^'^). Their weak limit exits and 
is exactly equal to Dh, because the mappings converge uniformly to h. 
Therefore, 

/ - = / ('|D/i^|2 + |m|2-2(D/i^,m)) 

^2 [ (\Dh\''-{Dh2,Dh)) 

= 2 f {Dh, Dh - Dh"^) 

We can now determine the number n = n,^ of congruent dyadic squares in 
Q,^, simply requiring that 

{diam ^ e for every i = 1, 2, . . . , n,^ 
||m--m||^.(U;;)^e-2- 

Fix such n = Hi, and abbreviate the notation for h^" to h'^. We obtain a 
homeomorphism 

oo 

hi:=h + Y^lh" -h]oeh + j^o{^) 

where we recall that [h^ — h]o stands for the zero extension oi — h to 
the entire domain $7. Clearly, hi is harmonic in each UJ,, v = 1,2,..., 
i = \,2, . . . ,ny and we have 

11^1 - h\\%'(i^) ^ sup{diamQ^ : = 1, 2, . . . , z = 1, . . . , n,,} < e 

oo oo 

(3.5) - ^ e + J^P^" - Dh\\^2(jj^^ ^ e + ^ e • 2"^^ = 2e 

For further considerations it will be convenient to number the squares QJ, 
and their preimages using only one index. These sets will be respectively 
denoted by Q° and U"^, a = 1,2,.... For the record, 

(3.6) diamQ°^e, a = 1,2,... 
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Finite energy case. Summing up the energy inequahties for the mappings 
/ij, : — )• Q1 we see that the total energy of hi does not exceed the energy 
of h. Even more, since h was assumed to be not harmonic, there is at least 
one region UJ, for which /i: UJ, — )• QJ, was not harmonic. Consequently, its 
harmonic replacement results in strictly smaller energy. Hence 

(3.7) £n[hi] < £n[h], so let 5 = \\Dh\\ ^2 (^q^ - \\Dhi\\ ^2 (^q) > 

Step 2. Denote by J-" = {Q" : a = 1,2, . . .} the family of all open squares 
QQ i~ Qa i~ Q* g^j.g Step 1 for the construction of the mapping 

/ii : O — )• f2*. Let V be the set of vertices of these squares. Whenever two 
squares Q°',Q^ ^ J-, a ^ f3, meet along their boundaries the intersection 
= dQ^ndQl^ is either a point in V or a closed interval with endpoints in 
V. Denote by J C {/"'^ : (3, a, /3 = 1, 2, . . . } the subfamily of all such 
intersections, excluding empty set and vertices. For each interval I"'^ G J7 
we shall construct a doubly convex lens-shaped region L'^'^ with I'^'^ as its 
axis of symmetry in the following way. Let -R be a number greater than 
the length of I""^ to be chosen later. There exist exactly two open disks of 
radius R for which /"''^ is a chord. Let L^'^ be their intersection. This is a 
symmetric doubly convex lens of curvature Thus L^'^ is bounded by two 
circular arcs 7"'^ = Q° n 9L^'^ and 7^^'" = Q'^ n dlJ^^ . As the curvature 
of the lens approaches zero the area of L^'^ tends to 0. This allows us to 
choose R depending on a and /? so that the lenses L"'^ = ^'^^ have the 
following property. 



The lenses L"'^ are disjoint because the opening angle of each lens is at 
most 7r/3 and their axes are either parallel or orthogonal. However, the 
closures of the lenses considered here may have a common point that lies in 
V. On each K"'^ we replace hi by the harmonic extension of its restriction 
to dK."'^ . Thus we obtain a homeomorphism h2'^ : K"^'^ L"''^ of class 

hi + £/oiK"'''^). By Proposition O the mappings h^''^ : IK"'^ ^ L"''' are 
diffeomorphisms. Finally, we define 



(3.8) 




where W 



■a,P — 



hi\Kin 



h2 = hi + - hi]. £ hi + s^oi^) = h + ^oi^) 

a,f3 



and observe that, from (j3.6p 
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Also, ()3.8p and Dirichlet's principle imply 

\ \Dh2 - Dhxf \ 2 (|D/i9p + ^ 4 V / \Dh 



2 

ll 



«,/3 Q,/3=r 

oo 2 
a,/3=l 

Thus 

II /i2 -/ii|L(/(f7) ^ e + 2e = 3e 

The boundary of K"'^ consists of two "^""-smooth arcs F"''' and F'^'" 
which share common endpoints, called the apices of W^'^ . These are preim- 
ages of 7°^'^ and 7^^'° under the mapping h\, respectively. Outside of the 
apices, the homeomorphism /12 : K"'/^ L"'^ is smooth with positive 
Jacobian. The smoothness is a classical result of Kellogg; a harmonic func- 
tion with 'if °^ -smooth values on a smooth part of the boundary is 'if ^ -smooth 
up to this part of the boundary [TBI Theorem 6.19]. The positivity of the Ja- 
cobian on such part of the boundary follows from the convexity of its image, 
see P p. 116]. 

In conclusion, /12 is locally bi-Lipschitz in 17 \ h^^iy). The exceptional 
set h~^{y) is discrete because V is. 
Finite energy case. By ()3.7p we have 

(3.9) \\Dh2\Wz(^n) ^ ll^/ii||^2(f,) ^ ||m||^2(f,) - 5 

Step 3. First we cover the set of vertices V by disks {D^ : f G V} centered 
at V with radii small enough so that 

(3.10) diamD^ ^ e, 

and {3B„ : u G V} is a disjoint collection of disks in 0*. Moreover, their 
preimages under /12 must satisfy 

(3.11) V / 1^2^ < 

Denote by Q.* = Q.* \ \J^(,y^v and = \ [j^^y /i^^(l^). Our focus for a 
while will be on one of the circular sides of a lens L°'^, say 

We truncate it near the endpoints by setting 7°'^ = O n 7"'^. Such trun- 
cated open arcs are mutually disjoint; even more, their closures are iso- 
lated continua in Q*. This means that there are disjoint neighborhoods of 
them. We are actually interested in a neighborhood of 7"'^ of the shape 
of a thin concavo-convex lens that we shall denote by L"'^. By definition, 
7"'^ C L"'^ C Q". The construction of such lens goes as follows. Let 
a and b denote the endpoints of 7""^, we assemble two circular arcs 7"'^ 
and 7"'^ with endpoints at a and b to form together with their endpoints a 
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concavo-convex Jordan curve. This Jordan curve constitutes the boundary 
of a circular lens L"'^. The term concavo-convex lens refers to the config- 
uration in which L"'^ lies in the concave side of the arc 7"''^ and convex 
side of 7"'^. It is clear that such lenses can be made arbitrarily thin so that 
L"''^ C 0* and the closures of L"'^ will still be isolated continua in Q*. From 
now on we fix the family {L'^''^ : a 7^ /3} of such concavo-convex lenses asso- 
ciated with the arcs 7"'''. We then look at their preimages U"''^ = ^(L"''^) 
and the "^""-smooth arcs T"^'^ = "^(7'^'^). The endpoints of T"^'^ lie in 
c?U"'^. Moreover, T"'^ splits U"''^ into two disjoint subdomains U"'^ and 
U"''^ such that U"'^ \ T"'^ = U"''' U U"'''. Here we have a homeomor- 
phism /i2 : U"''^ L"''^ which is '^°°-diffeomorphism on U"'^ and 
diffeomorphism on U_' . Therefore, for some positive number M^^p, we 
have pointwise inequlities \Dh2\ ^ M^^p and detL'/i2 ^ jj" — in both U"''' 

and U_' . Having established such a deformation of lenses and their preim- 
ages under /12 , we apply Corollary 15. 4[ We infer that there is also a constant 
M'^ p > with the following property: to every neighborhood of T"'^, say 

an open connected set Uo'^ C U"''^ that contains T"'^, there corresponds a 
"^""-diffeomorphism, denoted by h^: U"''^ L"'^, such that 



(3.12) 



'h^{z) = h2{z) for z E U"''' \ Uo ''^ 



\Dh3\ ^ M' and detD/is ^ -j-^ in U"'/^ 



We emphasize that ^ is independent of the neig hborhood Uo'^. We 
choose and fix Uq thin enough to satisfy 

• SUp^Q,/3|D/l2| ^ Ma^p 

• in the finite energy case, we also assume that |Uo'''| ^ [M^ 4~"~'^~^ 

Recall that 5 was defined by ()3.7p and later appeared in (j3.9p . This is 
certainly possible; for instance, take Uo''^ to be the preimage under /i2 of a 
sufficiently thin concavo-convex lens containing 7"'^. We call /13 : U"'^ 
L"'^ a smoothing of /12 : U"''^ L*^'^ associated with a given arc T"'^ = 
^2"H7°'^)- We now define a homeomorphism /13 : Q, — > Q.* by the rule 



h 



3 



{smoothing of /i2 in U°'^ 
h2 
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It belongs to /12 + =2^(0). Obviously /13 is a '^°°-diffeomorphism in We 
have for every z G 0, 



\h3iz) - h2iz)\ ^ 



diamL"'^ for z G U"''' 
otherwise 
^ diamO" ^ e 



see (j3.6p . Hence ||/i3 — /i2||<r(n) ^ £• As regards the energy of /is — /i2 we 
find that 

£n[h3-h2] = V / \Dh3-Dh2f 

.01 o^ ^J^IK'^I sup (1^31 + 1^21)2 

00 2 

These estimates sum up to 

\\h3 - h2\\^(n) ^ e + e = 2e 

Let us record for subsequent use the following estimate, obtained from (j3.1ip 
and (I3J3]1 . 



[ \Dhs\' + 2Y,([ 



|Z)/l3|2+ / ID/13P 



<l iDhsl' + iy il \Dhs-Dh2\'^+ \Dh2''^ 



3 — un2\ -r I l-L^ml 

) Jh2\30^) 

(3.14) ^ Yl i^ipf l^o'^l + + 2e' ^ Se^ 

Finite energy case. For the energy of /is, we observe that 

Il^/l3||i?2(f7) ^ ll^^3||^2(f^\uU?''') +Ell^^3||^2(u-.'5) 
^ P/l2||^2(j,) + J]|U^'^|l/2 SUp|D/l3| 



Thus 



(3-15) ll^^3||^2(f^) ^ ||-D/i||if2(j^) 



5_ 
2 
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Step 4. We have already upgraded the mapping /i to a homeomorphism 
/13 : — )• 0* such that /13 S /i + £/o{^) and 

/o -i^N 11^3 - Hj^{n) ^ 11^3 - ^2|U{n) + 11^2 - ^i||£/(n) + 11^1 - ^lU(n) 
(^•'^^ ^2e + 3e + 2e = 76 

Moreover, /i3 is a 'r°°-diffeomorphism on J7 \ U^gv h2^{B^). We now define 

a homeomorphism : 17 fi* by performing harmonic replacement of 
on each set h'^^{2S}y). This gives us a ^°°-diffeomorphism h^: h^^{20^) 
2B„, see Step 2 for details. For each z £ 

, , , , f2diamD^ if z G /i^^(2D^) 
I U otherwise 

Hence ||/i4 — /i3||<^(n) ^ 2e. Using ()3.14p we estimate the energy as follows. 

vev 



Thus, by (imi) 

(3.17) ||/i4-/i3lU(n) ^e + ^^e^6e 

Finite energy case. By virtue of Dirichlet's principle and (j3.15p we have 

(3.18) ll-D/i4||j?2(n) ^ ll-D/i3||^2(f^) ^ \\Dh\\^2(^^) - - 

Step 5. The final step consists of smoothing in a neighborhood of 
each smooth Jordan curves Cy = S/ig ^(2D^). We proceed in much the same 
way as in Step 3, but we appeal to Corollary 15.51 instead of Corollary 15.41 
By smoothing in a sufficiently thin neighborhood of each Cy we obtain a 

"^"^-diffeomorphism : Q Q*, /15 G /14 + ^o{^) such that 

(3.19) II/15 - /i4|U(n) ^ e 

We now recapitulate the estimates (|3.16p , (|3.17p and (|3.19p to obtain a 
diffeomorphism in Q 

H ■.= h^e /i + i/(S7) 

such that 

ll-f^ - ^IU(Q) ^ 11^5 - /i4|U(n) + 11^4 - h^\\si/(n) + 11^3 - ^lU(n) 

^ e + 6e + 7e = 14e 
which is as strong as (ii) in Theorem 11.21 

Finite energy case. To obtain the desired energy estimate £n[h^] ^ £n[h], 
we need to sharpen the energy part in (j3.19p . By narrowing further the 
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neighborhoods of C„ we can be make the energy £n[hc^ — /14] as small as we 
wish; for example to obtain 

\\Dh5 - -D/i4||^2(n) < - 
This is enough to conclude that 

because of (|3.18p . □ 

4. HOPF DIFFERENTIALS, THEOREM 11.31 

A quadratic differential on a domain in the complex plane C takes the 
form Q = F{z)dz (8) dz, where -F is a complex function on Q. Given a 
conformal change of the variable z, z = where ip: Q' ^ Q, the pull 

back 

ip^Q) = F{ip{0) dip dip = F{ip{0)'pHO «> 
defines a quadratic differential on Q'. It is plain that for a complex harmonic 
function /i : — t- C the associated Hopf differential 

Qh = hzhzdz (g) dz 

is holomorphic, meaning that 

(4.1) ^ {hj;} = 



dz 

Conversely, if a Hopf differential Qh = h^hz dz®dz is holomorphic for some 
■^^-mapping /i, then h is harmonic at the points where the Jacobian determi- 
nant J{z,h) := detDh = |/i2p — j/ifP / 0, see |T0', 10.5] and our Remark [Ol 
Here the assumption that J{z, h) ^ is critical. Let us illustrate it by the 
following. 

Example 4.1. Consider a mapping h G ^^'^(Co) defined on the punctured 
plane Co = C \ {0} by the rule 

for < l^l ^ 1 
+ I) for 1 ^ l^l < oo 



(4.2) h{z) -- 

Direct computation shows that 




and 



Thus 



hz{,z) 




for < |z| ^ 1 
for 1 ^ |z| < oo 

for < |2;| ^ 1 
for 1 ^ \z\ < oo 



I A n\ ^ dz (g) dz 
(4.3) Qh = ^ 
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It may be worth mentioning that the mapping h in (14. 2|) is the unique (up 
to rotation of z) minimizer of the Dirichlet energy 

g[H] = [ \DH\'^ 

over the annulus A = A{r,R) = {z: r < |2;| < R}, < r < I < R, 

subject to ah weak hmits of homeomorphisms H : A A{1,R^,), where 
i?* = ^ (iZ + ^) , see Note that the Hopf differential of (|4.3p is real along 
the boundary circles of A. The concentric circles are horizontal trajectories 
of Qh- In fact this is a general property of minimizers |21i Lemma 1.2.5]. 
The general pattern is that with the loss of injectivity comes the loss of the 
Lagrange-Euler equation for the extremal mapping. 

Properties of the function h with holomorphic Hopf differential Q = 
hzhz dz ® dz are of interest in the studies of harmonic mappings [HI 
[2T| [30| [3T] , minimal surfaces [51 [H^ and Teichmiiller theory [15] . In this sec- 
tion we prove Theorem 11.31 which imposes fairly minimal assumptions that 
imply harmonicity of y^^'^-solution to the equation ()4.ip . Some elements of 
the proof go back to [281 EH] . 

Proof of Theorem \1.3[ As a consequence of the Stoilow factorization theo- 
rem [H p. 56] the branch set of h is discrete, hence removable for continuous 

harmonic functions. Thus we assume that h: Cl Q* is a homeomorphism 
of Sobolev class such that 

(4.4) hzhz = F{z) is holomorphic in 

By virtue of Theorem ll.2| there exists a sequence of diffeomorphisms W : Q — 
Vl* converging c-uniformly and strongly in to h. Denote by 

(4.5) hih^ =: G J^Lm 

Thus F^ — )• F strongly in ^[^^(fi). Let us first dispose of an easy case. 
Case 0. The homogeneous equation F = 0. Since h^ are diffeomorphisms 
the Jacobian determinant J(z,h^) = — is either positive ev- 

erywhere in or negative everywhere in 0. Let us settle the case when 
J{z, h^) > for infinitely many indices j = 1, 2, ... . For such j we have 
\hi\ > I /ill, which yields |/i|p ^ \hihl\. Passing to the .if-'^-limit we obtain 

\h-z\^ ^ \hzh-z\ = \F{z)\=0. 

Thus h is holomorphic, by Weyl's lemma. Similarly, in case J{z, h^) < for 
infinitely many indices j = 1, 2, . . . , we find that h is antiholomorphic. 

Remark 4.2. We observe, based on the above arguments, that for this ho- 

1 2 

mogeneous equation hzhz = every solution h G (Q) obtained as the 
weak /^^'^-limit of homeomorphisms is either holomorphic or antiholomor- 
phic. The situation is dramatically different if hzhz ^ 0; some topological 
assumption on h are necessary, as illustrated in Example I4.1[ 
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Case 1. Nonhomogeneous equation F ^ 0. The function F, being holo- 
morphic, may vanish only at isolated points. Since isolated points are 
removable for bounded harmonic functions, it suffices to consider the set 
where F ^ 0. Proceeding further in this direction, we may and do assume 
that F{z) = 1 (by a conformal change of the z-variable) and /i is a W^"^- 
homeomorphism in the closure of the unit square Q = {3; + iy:0<x< 
1,0 < y < 1}. The problem now reduces to establishing that the equation 

(4.6) hzhl = 1 

implies A/i = 0. This will be proved indirectly by means of the energy- 
minimizing property 

(4.7) £Q[h]^£Q[H] 

where : Q — t- /i(Q) is any homeomorphism m. h + ^(Q); in particular, 
H = h on dQ. Indeed, if h were not harmonic, we would be able to decrease 
its energy by harmonic replacement (Propositions 12.11 and 12. 3p , contradict- 
ing (I32D. 

4.1. Proof of the inequality ()4.7p . With the aid of the approximation 
theorem we need only prove (j4.7p for mappings H £ h + £/o (Q) that are dif- 
feomorphisms on Q. From now on we assume that this is the case. Denote 
Q* = h{Q) = H{Q). We consider a sequence £ h + ^o(Q) of diffeo- 

morphisms : Q Q* converging in £/ (Q) to h. Moreover we may also 
assume that Dh^ — t- Dh almost everywhere in Q by passing to a subsequence 
if necessary. Now the sequence X"* • Q ~^ Q of self-homeomorphisms of the 
closed unit disk given hy X'' = o , where X'^ = id on dQ, is converging 

— 1 12 

uniformly on Q to x = H~ oh. It is important to observe that x ^ ('^) 
and x^ converges to x ii^ ^""^'^(Q') on any compactly contained subdomain 
Q' <s Q. Since and (x"')""'^ a-i'e diffeomorphisms on Q' and X''(Q')) respec- 
tively, the chain rule can be applied to the composition H = o (x"')"^- 
For w G X"'(Q') have 



dw ^ dw ^ dw 



dw ^ dw ^ dw 

where z = (x-^)~^(tt;). 

The partial derivatives of (X'' )^"'^ at ^ can be expressed in terms of xi 
and xi(^) by the rules 

djx'r' ^ xijz) 

dw J{z,x^) 
dw J{z,X^) 
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where the Jacobian determinant J{z,x^) is strictly positive. This yields 

dH ^ hixi - Hxi 
dw J{z,x') 

dw J{z,x^) 

We compute the energy integral of H over the set X"'(Q') by substitution 

w = x^{z), 



\hixi - hlxl? + \hixi - hixi? 



On the other hand, the energy of over the set Q' is 



Subtract these two integrals to obtain 

ViP + |/iiP)-|xiP-2Re 



hihlxixi 



dz 



(4.8) 



>4 



2\hihi\\xi\^-2Re 



hihlxixi 



dz 



\xi? - IxiP 



\hihU dz 



where we have introduced the notation 



1 



if hihi / 
otherwise. 



Note that Icr-'l = 1 and cr-' — t- 1 almost everywhere. 

Upon using Holder's inequality we continue the chain (j4.8p as follows. 

1 2 



(4.9) 



> 4 



/o' 



Xz - cr^Xz 



\hihi\dz 



Iq' J{z,hj)dz 
The denominator in ( 14. 9p is at most 1 because 

j{z,h^)dz = \x'm\ ^ 



Hhi\. 



1. 
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Therefore, 

Sq[H]-SQ\h^]^A 



hihil dz 



IM/iil dz. 



It is at this point that we can pass to the hmit as j — )■ cxd, to obtain 

1 2 



(4.10) 



£q[H] - £Q>[h] ^ 4 



\Xz 



Xz\ d2; 



Since Q' was an arbitrary compactly contained subdomain of Q, the esti- 
mate (|4.10p remains vahd with Q' replaced by Q. 

n 2 



(4.11) 



^4 



dx 

— dxdy 
dy 

^ dx{x,y) 





dy 



dy 



dx 



4 [ \xix, 1) - xix, 0)1 dx - 4 = 4- 4 = 
Jo 



as desired. 



□ 



Remark 4.3. When specialized to the case h G Theorem II .31 shows that 
h is harmonic outside of the zero set of its Jacobian. 

5. Auxiliary smoothing results 

Here we present some results concerning smoothing of piecewise differen- 
tiable planar homeomorphisms. They can be found in [25] in greater gener- 
ality, but since we require quantitative control of derivatives, a self-contained 
proof is in order. Here it is more convenient to use the operator norm of a 
matrix, denoted by ||-||. Note that ||j4|| ^ |j4| ^ 2||A|| for 2 x 2-matrices. 

Proposition 5.1. Let U C M? be a domain containing an open segment I 
with endpoints on the boundary dV which splits U into two subdomains Ui 



and U2 such that U \ / = Ui U U2 
homeomorphism with the following properties: 



Suppose that f : U U* 



C 



IS a 



the restriction of f to Vj is 'if °° -smooth, equals the 

1,2 the restriction of f 



i%) For 3 = 1,2,. 

identity on I; 
(a) There is a constant M > such that for j 
to U~ satisfies \\Df\\ ^ M and det Df ^ M'^ 

Then for any open set Uo with I C Uo C U there is a 'if°°-diffeomorphism 
g: U U* such that 

• g agrees with / on U \ Uo ( and also on I ); 

• \\Dg\\ ^ 20M and det Dg ^ (20Af)~i on U. 

Proof. Without loss of generality / C M = {{x,y): y = 0}. We write / 
in components as (n, v) where u and v are functions of x and y. Let us 
introduce a notation; given any "^""-smooth function /3: M — >• [0,oo), denote 
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^W) = ^ \y\ < /^{x)}- We can and do choose /? so that I C 

V{P) C Uo, and further scale it down until the following holds. 

I/5'(^)K4^ forallxGR; 
(5.1) |^;^| ^ _L_ mV{(3)\I, because v(x, 0) = 0; 

\ux ~ ^ Yo ™ ^(i^) \ because u{x,0) = x. 



As a consequence of (jn]) and (jS.ip . 

Since u is also M-Lipschitz by ([n]), the following double inequality holds in 
Vi(3) \ I. 

(5.3) ^ ^ - ^ M. 

Let us fix be a nondecreasing "^"^ function q: M — )■ M such that a{t) = 
for t ^ 1/3. Let a{t) = 1 for t ^ 2/3. Moreover, «'(*) ^ 4 for all t E M and 
a(oo) = 1, by convention. Now we introduce a modification of u on U by 
setting 

u := a{t)u + (1 — a[t))x where ^ = < 

I oo otherwise 

Note that u = u outside of V{/3). In V{I3) \ / we compute the derivatives as 
follows. 

= -t'^a'{t)f3'{x)'^^^-^ + a{t)ux + 1 - a{t) 

(5.4) 

Zi X 

Uj, = ta'{t) h a{t)uy 

Since u is M-Lipschitz by (jn]), we have |m — x| ^ -M|?/|. From this, (|5.ip 
and (15.41) we obtain 



(5.5) — ^ — , and \uy\ ^ 5M, 



8 12 

^ Ur ^ 

10 10 

which combined with (j5.2p yields 

8 1 _ 5M _ 3 
(5.6) n^t;, - ^ ^0 2M " SOM^ " lOM' 

Next we modify v on U. Specifically, 

V := a(sjw + (1 — CK(sjj77T7 where s = < ''^^ 

2M I oo otherwise 

Note that v = v outside of V{(5/2>), and on the set y(/3/3) we already have 
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Computations similar to ()5.4|] yield (on the set V{l3/3) \I) 
Vx = --a {s)s — — + a[s)vx] 

sa'is) ( y \ , , 1 - a{s) 

Vy = — iv —] +a(s)vy H —. 

y y \ 2MJ ^ ' ^ 2M 

Straightforward estimates based on (|5.ip . ()5.2p and (|5.3p comply 

4M 1 3M 

(5,8) 'I'' * — Bom? * — ■ 

It remains to check that the mapping g := {u,v), which agrees with / 
outside of V{/3), satisfies all the requirements. As regards '^""-smoothness 
we need only check it on V{f3/9). But in this neighborhood of / we have 
a linear mapping, g{x,y) = {x,^), so '^""-smooth. By virtue of (jS.Sp 
and ()5.8p we have \\Dg\\ ^ 20M. The desired lower bound for detDg follows 
from ()5.6p and (jS.Sp . Consequently, g is a local diffeomorphism, and since it 
agrees with / on dV{l3), it is in fact a diffeomorphism, by a topological result: 
a local homeomorphism which shares boundary values with a homeomorphism 
is injective \25{ Lemma 8.2]. □ 

We also need a polar version of Proposition 15.11 

Corollary 5.2. Let U C 6e a domain containing a circle T. Suppose 
that / : U ^ U* C is a homeomorphism with the following properties: 
(i) The restriction of f to T is the identity mapping; 

(a) There is a constant M > such that the restriction of f to either 
component o/U\T is '^'^ -smooth with \\Df\\ ^ M anddetDf ^ M^^. 

Then for any open set W with T C W C V there is a -diffeomorphism 
g : U — )• U* such that 

• g agrees with f onV\W and on T; 

• \\Dg\\ ^ 80M and detDg ^ (80M)^i on V. 

Proof. It is convenient to identify with C. Without loss of generality 
T = {z £ C: \z\ = 1}. Let ip{C) = exp(i(^). The mapping F = ip~^ o f o ip is 
well-defined in some open horizontal strip Sh = {z £ C: \lmz\ < e} which 
we choose thin enough so that 1^(8^) C W and \ip'\'^ < e^*^ ^ 2. Note that F 
is 27r-periodic and satisfies 

\\DF\\ ^ 2M and det DF ^ (2M)~^ 

The proof of Proposition 15.11 applies to F with no changes other than one 
simplification: /3 > is now a small positive constant rather than a function. 
Thus we obtain a diffeomorphism G which agrees with F on MU (5 \ V{/3)) 
and satisfies \\DG\\ ^ 40 and det DG ^ (40M)-^ Since F was 27r-periodic, 
so is G. Thus, g := ip o G o is the desired diffeomorphism. □ 
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Our applications require slightly more general versions of Proposition 15.11 
and Corollary l5.2l where the separating curve is allowed to have other shapes 
and / is not required to agree with the identity on the curve. 

Definition 5.3. A parametric curve T: (0, 1) — t- is regular if T extends 
to a bigger interval (a, b) D [0, 1] so that the extended mapping is a 
diffeomorphism onto its image. 

Note that a regular curve F has well-defined endpoints r(0) and r(l). 
Also, r extends to an injective "^°°-mapping <1>: (0, 1) x (—1, 1) — )• M? such 
that ll-D*!?!! and ||(Z)<I>)~^|| are bounded. This follows from the existence of 
a tubular neighborhood of the image of T [23j Theorem 4.26]. 

Corollaries 15.41 and 15. 5| given below, generalize Proposition 15.11 and Corol- 
lary 15.21 respectively. 

Corollary 5.4. Let U C M? be a domain containing the image of a regular 
arc T with endpoints on the boundary dV which divides U into two subdo- 
mains Ui and U2 such that U \ P = Ui U U2 . Suppose /: U IF C is 
a homeomorphism such that f oT is also regular and the restriction of f to 
each Uj is '^'^ -smooth and satisfies 

\Df{z)\ ^ M, deiDfiz) ^ ^ forz^Vi 

where M is a positive constant. Then there is a constant M' > such that to 
every open set U' C U with P C U' there corresponds a -diffeomorphism 
g : U U* with the following properties 

• g{z) = f{z) for z € U \ U' ( and also onV) 

• \Dg{z)\ ^ M' and deiDg{z) ^ onV. 

Proof. Let Q = (0, 1) x (-1, 1). Let <i> and ^ be the extensions of P and foT 
to Q as in Definition [Ol There is a domain U such that (0, 1) x {0} C_U C Q, 
$(U) (s U', and the composition F := o / o $ is defined in U. Note 
that F = id on (0, 1) x {0}. We apply Proposition 15. II (with U in place of U 
and with F in place of /) and obtain a '(^°°-diffeomorphism G: U — )• F{V). 
Finally, replace F within U with the diffeomorphism g = o G o . □ 

Corollary 5.5. Let U C 6e a domain containing the image of a 'io°°- 
smooth Jordan curve P which divides U into two subdomains Ui and U2 such 
that U \ P = Ui U U2 . Suppose / : U — )• U* C is a homeomorphism such 
that the restriction of f to each Uj is 'io°° -smooth and satisfies 

\Df{z)\ ^ M, deiDf{z) forzGVi 

where M is a positive constant. Then there is a constant M' > such that to 
every open set U' C U with P C U' there corresponds a -diffeomorphism 
g : U U* with the following properties 

• g{z) = f{z) for z G U \ U' ( and also onV) 
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• \Dg{z)\ ^ M' and det Dg{z) ^ ^ onU. 
Proof. The proof of Corollary 15.41 is easily adapted to this case. □ 

6. Concluding remarks 

One may wonder whether the proof of Theorem 11.21 can be extended to 
the spaces W^'"^, 1 < p < oo, by means of the p-harmonic replacement in 
place of Proposition [231 Indeed, p-harmonic mappings are '"-smooth |36j. 
However, the injectivity of p-harmonic replacement of a homeomorphism is 
unclear. 

Question 6.1. Is there a version of the Rado-Kneser-Choquet theorem for 
p-harmonic mappings? That is, does the p-harmonic extension of a homeo- 
morphism onto a convex Jordan curve enjoy the injectivity property? 

An attempt to extend Theorem 11.21 to higher dimensions faces another 
obstacle: the Rado-Kneser-Choquet theorem fails in dimensions n ^ 3 as 
was proved by Laugesen [22] . 

References 

1. K. Astala, T. Iwaniec, and G. Martin, Elliptic partial differential equations and qua- 
siconformal mappings in the plane, Princeton University Press, Princeton, NJ, 2009. 

2. J. M. Ball, Global mvertibility of Sobolev functions and the interpenetration of matter, 
Proc. Roy. Soc. Edinburgh Sect. A, 88 (1981), 315-328. 

3. J. M. Ball, Discontinuous equilibrium solutions and cavitation in nonlinear elasticity, 
Philos. Trans. R. Soc. Lond. A 306 (1982) 557-611. 

4. J. M. Ball, Singularities and computation of minimizers for variational problems. 
Foundations of computational mathematics (Oxford, 1999), 1-20, London Math. Soc. 
Lecture Note Ser., 284, Cambridge Univ. Press, Cambridge, 2001. 

5. P. Bauman, D. Phillips, and N. Owen, Maximum principles and a priori estimates for 
an incompressible material in nonlinear elasticity, Comm. Partial Differential Equa- 
tions 17 (1992), no. 7-8, 1185-1212. 

6. J. C. Bellido and C. Mora-Corral, Approximation of Holder continuous homeomor- 
phisms by piecewise affine homeomorphisms, Houston J. Math., to appear. 

7. S. Conti and C. De Lellis, Some remarks on the theory of elasticity for compressible 
Neohookean materials, Ann. Sc. Norm. Super. Pisa CI. Sci. (5) 2 (2003) 521-549. 

8. U. Dierkes, S. Hildebrandt, A. Kiister, and O. Wohlrab, Minimal surfaces. I. Boundary 
value problems. Springer- Verlag, Berlin, 1992. 

9. P. Duren, Harmonic mappings in the plane, Cambridge Tracts in Mathematics, 156. 
Cambridge University Press, Cambridge, 2004. 

10. J. Eells and L. Lemaire, A report on harmonic maps. Bull. London Math. Soc. 10 
(1978), no. 1, 1-68. 

11. J. Eells and L. Lemaire, Selected topics in harmonic maps, CBMS Regional Conference 
Series in Mathematics, 50. American Mathematical Society, Providence, RI, 1983. 

12. J. Eells and L. Lemaire, Another report on harmonic maps. Bull. London Math. Soc. 
20 (1988), no. 5, 385-524. 

13. L. C. Evans, Quasiconvexity and partial regularity in the calculus of variations. Arch. 
Rational Mech. Anal. 95 (1986), no. 3, 227-252. 

14. L Fonseca and W. Gangbo, Local invertibility of Sobolev functions, SIAM J. Math. 
Anal. 26 (1995), no. 2, 280-304. 



HOPF DIFFERENTIALS AND SMOOTHING SOBOLEV HOMEOMORPHISMS 21 



15. L. Gardiner and N. Lakic, Quastconformal TeichmiiUer theory, American Mathemat- 
ical Society, Providence, RI, 2000. 

16. D. Gilbarg and N. S. Trudinger, Elliptic partial differential equations of second order, 
Springer- Verlag, Berlin-New York, 1977. 

17. F. Hcloin, Homeomorphismes quasi conformes entre surfaces riemanniennes, C. R. 
Acad. Sci. Paris Sr. 1 Math. 307 (1988), no. 13, 725-730. 

18. S. Hcncl and P. Koskela, Regularity of the inverse of a planar Sobolev homeomorphism. 
Arch. Ration. Mech. Anal. 180 (2006), no. 1, 75-95. 

19. H. Hopf, Differential geometry in the large. Notes taken by Peter Lax and John Gray. 
With a preface by S. S. Chern. Lecture Notes in Mathematics, 1000. Springer- Verlag, 
Berlin, 1983. 

20. J. Jost, A note on harmonic maps between surfaces, Ann. Inst. H. Poincare Anal. Non 
Lineaire 2 (1985), no. 6, 397-405. 

21. J. Jost, Two-dimensional geometric variational problems, John Wiley & Sons, Ltd., 
Chichester, 1991. 

22. R. S. Laugesen, Injectivity can fail for higher- dimensional harmonic extensions Com- 
plex Variables Theory Appl. 28 (1996), no. 4, 357-369. 

23. S. Montiel and A. Ros, Curves and surfaces, American Mathematical Society, Provi- 
dence, RI, 2005. 

24. C. Mora-Corral, Approximation by piecewise affine homeomorphisms of Sobolev home- 
omorphisms that are smooth outside a point, Houston J. Math. 35 (2009), no. 2, 
515-539. 

25. J. Munkres, Obstructions to the smoothing of piecewise-differentiable homeomor- 
phisms, Ann. of Math. (2) 72 (1960), 521-554. 

26. S. Miiller, S. J. Spector, and Q. Tang, Invertibility and a topological property of Sobolev 
maps, SIAM J. Math. Anal. 27 (1996), no. 4, 959-976. 

27. Ch. Pommerenke, Boundary behaviour of conformal maps. Springer- Verlag, Berlin, 
1992. 

28. E. Reich and K. Strcbcl, On quasiconformal mappings which keep the boundary points 
fixed. Trans. Amcr. Math. Soc. 138 (1969), 211-222. 

29. E. Reich and H. R. Walczak, On the behavior of quasiconformal mappings at a point. 
Trans. Amer. Math. Soc. 117 (1965), 338-351. 

30. R. Schoen, Analytic aspects of the harmonic map problem, in "Seminar on nonhnear 
partial differential equations" (Berkeley, Calif., 1983), 321-358, Math. Sci. Res. Inst. 
Publ., 2, Springer, New York, 1984. 

31. H. C. J. Sealey, Harmonic diffeomorphisms of surfaces, in "Harmonic Maps", Lecture 
Notes in Mathematics vol. 949, 140-145, Springer, New York, 1982. 

32. G. A. Seregin and T. N. Shilkin, Some remarks on the mollification of piecewise-linear 
homeomorphisms. J. Math. Sci. (New York) 87 (1997), no. 2, 3428-3433. 

33. J. Sivaloganathan and S. J. Spector, Necessary conditions for a minimum at a radial 
cavitating singularity in nonlinear elasticity, Ann. Inst. H. Poincare Anal. Non Lineaire 
25 (2008), no. 1, 201-213. 

34. M. Struwe, Plateau's problem and the calculus of variations, Princeton University 
Press, Princeton, NJ, 1988. 

35. V. Sverak, Regularity properties of deformations with finite energy, Arch. Rational 
Mech. Anal. 100 (1988), no. 2, 105-127. 

36. N. N. Ural'ceva, Degenerate quasilinear elliptic systems. Zap. Naucn. Sem. Leningrad. 
Otdel. Mat. Inst. Steklov. (LOMI) 7 (1968) 184-222. 



22 TADEUSZ IWANIEC, LEONID V. KOVALEV, AND JANI ONNINEN 

Department of Mathematics, Syracuse University, Syracuse, NY 13244, USA 
AND Department of Mathematics and Statistics, University of Helsinki, Fin- 
land 

E-mail address: tiwaniecSsyr . edu 

Department of Mathematics, Syracuse University, Syracuse, NY 13244, USA 
E-mail address: IvkovaleSsyr . edu 

Department of Mathematics, Syracuse University, Syracuse, NY 13244, USA 
E-mail address: jkoniiine@syr.edu 



